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ABSTRACT 

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational 
implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical 
results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic 
in time and find that this also applies to narrow jets with half opening angle of 0.05 radians. As a result, 
afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways 
expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the 
jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to 
take the observer angle into account can lead to an overestimation of the jet energy by up to a factor 4. This 
weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late time 
radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is 
approximately on-axis and where differences of a few in flux level between the model and the simulation are 
acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high 
observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant 
late time features such as the rise of the counterjet, detailed jet simulations remain indispensable. 



1. INTRODUCTION 

Gamma-ray burst (GRB) afterglows are often theorized to 
result from synchro tron radiation from a d ecelerating rela- 
tiv istic blast wave (Meszaios & Reei ll997i see also reviews 
by Meszaios 20 061 |Piranll2005l and references therein). Over 
the past few years, different groups have performed two- 
dimensional (2D) relativistic hydrodynamics (RHD) GRB 
afterglow jet simulations at very high resolution, coupled 
to sy n chrotron radiation calcu l ations {Zhang & MacFadven 
2009; Van Eert en et al.1 l20lol l20Tlt Iwvgoda et al.1 1201 U 
van Eerten & MacFadven 2011). The two main reasons 
to employ RHD simulations are that both jet decollima- 
tion and deceleration are very difficult to capture in de- 
tail in analytical models. While self-similar solutions ex- 
ist for both the ultra-relativistic Blandford-McKee (BM) 
phase (Bland ford & McKedfl976l assuming a conic section 
of this spherical solution can be applied before the jet edges 
come into causal c ontact) and the late nonrel ativistic Sedov- 
Taylor (ST) phase (iTavloril 19501 ISedovl[l959l) . no fully self- 
consistent model exists connecting both regimes. The recent 
numerical studies cited above have begun to systematically 
explore the transrelativistic behavior of afterglow jets. In this 
paper we address the following questions regarding jet dy- 
namics and the shape of light curves calculated directly from 
recent high resolution numerical simulations: 



1. Recent numerical st udies dZha ng & MacFadven 2009; 
IWygoda etaLl l20TTb agree on the absence of expo- 
nential sidew ays expansion of the jet (as predicted by 
Rhoads 1999 to occur in the relativistic regime) for jets 
with a reasonably wide half opening angle of 0.2 ra- 
dians. Is this type of expansion significant for lateral 
expansion of narrower jets? 

2. Is th e logarithmic spreading (Zhang & MacFa dvenl 
12009!) common to both narrow and wide jets and what 
does that imply for transition to the non-relativistic 
regime? 



3. What do the answers to question 1 and 2 imply for the 
shape of the light curve as set by the competing effects 
of jet spreading and the jet edges becoming visible due 
to a decrease in relativistic beaming of the radiation? 

4. Will the observer angle still play a role for the observed 
jet break of narrow jets, as it does for wide jets? 

5. What are the implications for late time radio observa- 
tions, especially with respect to calorimetric estimates 
of the jet energy? 

Related to question 4 is the issue of the importance of the 
observer angle for the energy estimate of the GRB from 
the jet break. It has bee n claimed in observational studies 
dCeiikoet all 120 101 120111) that the magnetar model for GRB 
progenitors is severely challenged by recent observations (e.g. 
GRB 090926A). In view of this, we will also return to the 
question 

6. Is the magnetar model still challenged when the ob- 
server angle is taken into account in determining the 
jet energy from the jet break? 



that we briefly discussed in I Van Eerten et all (120101) . 

The rest of this paper is organized as follows. After a brief 
interlude establishing some notation and relevant timescales, 
in section [3] we discuss the dynamics of afterglow jets (i.e. 
questions 1, 2). In section|4]we discuss the observational im- 
plications (questions 3-6). We summarize our results in sec- 
tion 

2. TIMESCALES 

Analytical studies have argued for a number of timescales 
of special relevance. Some of these timescales will feature in 
the discussion below and we therefore define them here. 

The first of these is the point in time that approximately 
marks the onset of lateral spreading of the jet. During the 
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BM phase, and before any significant lateral spreading, the 
radial outflow decelerates according to 7 = 7o(?/*o)~ 3 > with 
7 the Lorentz factor of the fluid directly behind the shock. 
Here, 70 and ?o are determined completely by the explosion 
energy and circumburst medium density, the physical param- 
eters that set the scale of a BM point explosion. When the 
flow has decelerated to the point where 7 ~ 1/9q, where #0 
the original jet half opening angle in radians, sideways expan- 
sion is expected to take off. Two motivations for this are that 
this time tg marks approximat ely the time wh en the edges of 
the jet achieve causal contact (Shapiro 1979), as well as the 
point where the widening of the jet becomes comparable to 
the original jet opening angle (under the assu mption that th e 
relativistic jet spreads with the speed of sound, Rh oadsll999l) . 
Using conservation of energy in the expanding blast wave (eq. 
43, BM), we find 



tg =235(£,, o /lO 53 «o) 1/3 (0o/O.l) 2/3 days. 



(1) 



Here Ej S0 denotes the isotropic equivalent energy of the jet in 
ergs and «o the circumburst medium number density in cm -3 . 
Aside from the constant light travel time from the origin of 
the explosion to the observer, the observed time fy® differs 
significantly from the time in the progenitor frame because 
the jet almost keeps up with its own radiation. For the front of 
the jet at radius R, this transformation implies 



l =tg-R(tg)/c = 



tg 
I67 2 



(2) 



where we have used R = ct(l - I/I67 2 ) (eq 26, BM). Here c 
denotes the speed of light and 7 refers to the fluid Lorentz 

factor directly behind the shock fr ont. 

The outflow has been argued dLivio & Wa xman 2000) to 
become roughly spherical by the time t Sj ® (in the observer 
frame), expressed as 



t,.m = 230 



10 53 



1/3 



_1 
0.1 



V 3 



days, 



(3) 



in Wvgod a et al.l (1201 11) . Here E- lso denotes the isotropic 
equivalent energy of the jet in ergs and «o the circumburst 
medium number density in cm" 3 . If f 5) ® marks the point where 
the jet becomes quasi-spherical, a direct implication is that 
beyond this point in time the original orientation of the jet 
with respect to the observer has become irrelevant. Note that 
f.s.® ~ tg (the difference between 235 and 230 being a differ- 
ence in round-off), even though the two timescales are ex- 
pressed in different frames. The implication is therefore that 
the difference between frames has become neglible by this 
point in time. This is because two assumptions have been 
made in the derivation of t Si ®. it is assumed that lateral spread- 
ing is a very fast process that takes approximately tg ~ R(tg) / c 
to complete, and that during the spreading phase further in- 
crease in radius compared to R(tg) is neglible. Later we will 
see that the second of these assumptions leads to a grossly 
overestimated t , .0 compared to what fast spreading would im- 
ply, and that as a result t s $ in practice becomes comparable 
to the time when the onset of the nonrelativistic phase is ob- 
served. 

After this time f^®, the jet is expected (Piran 2005) to fur- 
ther decelerate and become nonrelativistic at time 



f NR =1100 



10 53 n / 



1/3 



days. 



(4) 
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FIG. 1 . — A snapshot of the comoving number density at t = 0.95 X (nr ~ 
9 X 10 7 s for the 6» = 0.05 rad, Ej = 1 .25 X 10 50 erg simulation (left) and the 
00 = 0.2 rad, Ej ■ = 2 X 10 5 ' erg simulation (right). Both simulations have the 
same £j so . The distances are in lightseconds. The outer radius along the axis 
of the narrow jet is 1 .9 X 10 18 cm and 2.55 X 10 18 cm for the wide jet, while 
c x t = 2.7 x 10 18 cm. 

It is worth noting that t s ,® depends on the original jet opening 
angle, whereas ?nr does not. We also n ote that the numerical 
factor in the equation for t s ,® given in Zhang & Ma cFadvenl 
(2009) is approximately tw ice that in [W vgoda et al. | (|201 ll) 
and the numerical factor in Zhang & MacFadven (2009) for 
?nr is 970 rather than 1100 days. We follow Wvgo da et al.l 
d201 ll) in order to compare directly to their work, but these 
differences do not alter our conclusions. Because at ?n R the jet 
has become nonrelativistic, it is assumed that c?nr 3> R(tNtd 
and that as a result ?nr m ?nr,©- 

3. JET DYNAMICS 

We have used the RHD code RAM (IZhang & M acFadven 
2006) to run a number of high-resolution jet simulations start- 
ing from the BM solution. A series of jets with 6*o = 0.05 rad 
and a wide jet with 9q = 0.2 rad, starting from fluid Lorentz 
factors 7 = 40 and 7 = 25 at the shock front Q and both with 
Eiso = 10 53 have been calculated as well as a 8q = 0.05 rad 
jet with an energy in both jets of Ej = 2 x 10 51 erg (£,,„ = 
2Ej/6q = 1.6 x 10 54 erg). In all cases, rip = 1 . The simula- 
tion re solution in all cases is similar to Zhang & MacFadven 
(2009) and t he same adaptive mesh refineme nt strategies as 
described in van Ee rten & MacFadvenl (1201 ll) have been ap- 
plied. 

Our simulations use an EOS which smoothly interpolates 
betwe en relativistic (T a d = 4/3) and non-relativistic (T a d = 
5 /3) dMignone et al.ll2005l) . The effect of the EOS on the dy- 
n amics and the shape of th e light curve are discussed in detail 
in lVan Eerten et akl (12010). There it was found (for spherical 
outflow) that keeping T a d = 4/3 led to a 9 percent smaller ra- 
dius at late time compared to a changing Y a d- The interpolat- 
ing EOS will lead to steeper light curves in the transrelativistic 
phase. 

1 Simulations with initial 7 > 1 /9q are expected to lead to nearly identi- 
cal fluid evolution (and light curves), as lateral spreading has not yet begun 
and the radial outflow is still expected to follow BM as confirmed by our 
simulations starting with different Lorentz factors 
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FIG. 2. — Top: Outer angle of three regions containing fixed fractions of 
total jet energy for the 9q = 0.05 rad, Ej = 1.25 X 10 50 erg jet. The vertical 
grey lines denote t(-y = 5), r(7 = 2) and ?nr from left to right. The dash-dotted 
black curve denotes exponential expansion using 9 = 8q exp[c(f — foVsNTL 
the dash-dotted blue curve the modified Rhoads model fit to f(7 = 5) and the 
dash-dotted red curve a logarithmic expansion curve fit to the 95 % curve 
between the onset of spreading and ?nr. Bottom: same as top, only rescaled 
using 0q — 0.05 and fg = f(7 = 1 /9) ~ 148 days, for direct comparison with 
Fig. 1 of Wygoda et al. 12011). The single vertical grey line denotes the 
scaled ?nr. 



Fig. [T]compares between narrow and wide jet at 0.95?nr 
10 3 days. The narrow jet has Ej = 1.25 x 10 50 erg, the wide 
jet Ej = 2x 10 51 erg (matching Fig. 2, lWvgoda et all (12011b ). 
The different Ej values lead to different rad ii around ?nr- 
When we compare our Fig. Q] to Fig. 2 in Wygo da et al.l 
(l20Tlh . we find excellent agr eement in the width a nd over- 
all radial extent of the jets in Wy goda et al.l d201 ll) and our 
simulations. 

We caution against interpreting the differences in radius for 
the wide and narrow jet visible in Fig. Q]as proof of fast early 
time spreading and deceleration of the narrow jet relative to 
the wide jet. Noticeable differences in radius between narrow 
and wide jet occur after 7 ~ 1 by definition and could there- 
fore be attributed to Ej rather than early time spreading or 
(9(). At this stage £i S0 no longer applies due to causal contact 
across all angles of the jet and neither does the assumption of 
ST self-similarity as long as the jet is not spherical and an ad- 
ditional length scale is introduced by the current width of the 
jet. 

We now examine the angular distribution of jet energy for 
the narrow jet with Ej = 1.25 x 10 50 erg. In Fig. [2] we have 
plotted the evolution of the boundaries marking the regions 



containing different fractions of the total jet energy. Before 
discussing the behavior of the 95 % boundary we emphasize 
that its interpretation as 'edge' of the jet will become ambigu- 
ous if the front of the jet is not homogeneous (i.e. when en- 
ergy and Lorentz factor depend on angle with respect to the jet 
axis). Thi s inhomogeneity is a key finding of previous numer- 
ical work dZhang & M acFadv en|2009tlVan Eerten et al.l201 ll 
iGranot et al .112001 I lGranoul2007l) . 

Taking the 95 % boundary at face value and comparing it 
with exponential expansion, we see from Fig. [2] that both 
curves start to deviate before 7^5 (for exponential expan- 
sion, t = frjexp[c(f - to) I Isnt], where to denotes the time when 
7 = 1/6*0, c the speed of light and Isnt ~ 2.7 x 10 17 cm the 
Sedov length for 60 = 0.05 rad., see Zhang & Mac Fadyenl 
(2009)). The 95 % opening angle at 7 = 5 is ~ 76* ~ 0.35 
rad is large with respect to &o, but still only a small fraction of 
the final opening angle, meaning that most lateral expansion 
will not take place in the strongly relativistic regime. 

The exponential curve can be improved upon by not tak- 
ing the expone ntial limit wh en implementing the dynamical 
equations from Rhoads 1999. We obtai n the dash-dott e d blue 
curve shown in Fig. [2] when following Wygoda et alJ (1201 1 ) 
by substituting d9(r)/r oc 1 /jr for eq. (3) in Rhoads ( 1999) 
and including an additional scaling factor determined by fit- 
ting to the simulation curve. However, logarithmic expansion 
provides a better fit, not only beyond 7 ~ 2 or 5, but for the 
entire region up to ?nr, even for narrow jets. The values for 
this fit function have been determined from 8 at ?nr and 8 at 
the onset of the sideways movement of the 95 % boundary 
(f ~ 150 days). Modified Rhoads ends in a logarithmic in- 
crease but is too constrained to get the slope right. 

Finally, Fig. [2] shows that the 95 % boundary lies at 0.7 
rad for the narrow jet at ?nr, but an extended phase of lateral 
expansion still follows during which the increase of the curve 
even temporarily becomes steeper than at earlier times. The 
non-relativistic dynamics depend on Ej rather than Ei SO and 
Ej jet rather than 6*o determines the late-time sideways expan- 
sion. In the bottom plot the angles are rescaled as fraction of 
the original opening angle and the time as fract ion of the time 
when 7 = 1/6*0, f° r comparison with Fig. 1. in Wygoda et al. 
d201 ll) . Note especially that the 90 and 50 % boundaries show 
no lateral spreading at 7 = 1/6*0 

However, as stated above, arguments invoking the 95 % 
boundary become less meaningfull if the shock front is in- 
homogeneous. Fig. [3] shows that at t = 0.95?nr (i-e. th e same 
time as in our Fig. Q]and Fig. 2 of lWygoda et al.1 ([201 1)) there 
is still an order of magnitude difference in energy between the 
50 % and the 95 % regions, directly confirming that the shock 
fronts remain inhomogeneous even on very long timescales. 
On the other hand, the 95 % boundaries for the wide and nar- 
row jets lie very closely together at this point in time, sug- 
gesting that at least superficially the jet retains less memory 
of 6*o and Ej is becoming the dominant factor determining jet 
dynamics. 

We conclude that even for narrow jets no phase of expo- 
nential spreading occurs that leads to approximate spheric- 
ity over only a logarithm ic increase in radius, as argued by 
iLivio & Waxmanl (yOOO) (note that t s <C ?nr): Exponential 
expansion does not fit the light curve down to 7 ~ 5 for 
95% and the 50% boundary only starts spreading in the tran- 
srelativistic regime. Instead, the expansion accross the en- 
tire transrelativistic r e gime is logarithmic, as described in 
Zhang & MacFadyen] d2009l) (and as analytically expected, 
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FIG. 3. — Angular distribution of energy at t = 0.95/nr for both wide (0.2 
rad, solid black curve) and narrow (0.05 rad, dash-dotted black curve) jet. 
The two vertical green lines on the left denote the outer angle of the region 
containing 50 % of the jet energy (again, solid for 0.2 rad, dash-dotted for 
0.05 rad). The two red lines on the right denoting the outer angle of the 
region containing 95 % are nearly indistinguishable. 
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FIG. 4. — X-ray afterglow light curves showing the jet break for the nar- 
row 0.05 jet with Ej = 2 X 10 52 erg typical to Swift light curves, for ob- 
servers at different angles. The calculation uses z = 0.5, luminosity dis- 
tance di = 8.7 X 10 27 cm, s ynchrotron slope p = 2.5 and the method from 
IZhang & MacFadven 12009). Vertical grey lines denote jet break fit results 
for Bbs = 0, O . 

see the appendix of IVan Eerten et ai]l2010l) . We find that the 
jets remain strongly inhomogeneous and highly nonspherical 
at late times. 

4. OBSERVATIONAL IMPLICATIONS 

4. 1 . jet break and observer angle 

Having established that even for narrow jets the jet dynam- 
ics are not dominated by exponential expansion, it follows that 
the observed jet break is also likely not dominated by lateral 
expansion of the flow. Even for narrow jets, features asso- 
ciated missing flux due to the jet edges becoming visible are 
still expected to play a role in, if not dominate, the shape of the 
ligh t curve. One such feature is the role of the observer an- 
gle dV an Eerten et al.ll201 U [van Eerten & MacFadven] 1201 It 
Van Eerten et al. 201 la). Fig. g]shows X -ray light curves (1.5 
KeV) for a narrow jet with Ej ■ = 2 x 10 52 erg seen from dif- 
ferent observer angles. Given that the observed jet break time 
tj was found to be set by the edge furthest from the observer, 



IVan Eerten et al.1 (120101) suggest 

= 3.5(1 J "o ( Q2 ) da y s > ( 5 ) 

where z denotes redshift. In practice jet spreading, radial 
fluid structure and arrival time effects render the light curve 
more complex than mere power laws. Smooth power law fits 
assuming a fractional error of 10 pet on the datasets from 
Fig. 2] lead to a jet break time 6.0 x 10 4 ± 4 x 10 3 s on- 
axis (x 2 /dof = 0.25) and 2.22 x 10 5 ± 6 x 10 3 s on-edge 
(X 2 1 do f = 0.26). 

Although the observer angle effect may therefore not be as 
severe as implied by eq. [5] this does confirm that the error 
due to observer angle in energy estimates based on jet open- 
in g angles is a genera l issue that not only applies to wide jets. 
In lCenkoet ail (1201 ll) GRB 090926A, with prompt energy re- 
lease E 1 = 1.4 x 10 52 erg, is found to exhibit the most severe 
challenge to the magnetar energy limit of 3 x 10 52 erg. To 
illustrate: downscaling the break time by a factor 3.7, down- 
scales the jet opening angle by a factor 1.6 and E y by a factor 
2.7, leading to £ 7 ~ 5.2 x 10 51 erg. Although we therefore 
conclude that the challenge posed to the magnetar model by 
events like GRB 090926A might be weaker than previously 
reported, we emphasize that we draw no conclusions regard- 
ing the valid ity of the magnetar model and that the significant 
finding from lCenko et al.1 (1201 ll) that a class of very energetic 
(relative to earlier observations) GRB events is emerging re- 
mains unaltered. 

If jet breaks for narrow and wide jets exhibit the same 
general shape, earlier conclusions regarding the effect of jet 
spreading for non-narrow jets are general as well. Specifi- 
cally, this implies that, while subdominant, jet spreading does 
affect the pos t-break slope. Th i s effe ct has been demon strated 
explic itly in IVan Eerten et al.l d201 ll) and earlier in iGranotl 
d2007l) . We note that jet expansion (even if logarithmic) can 
not be fully ignored when making quantitative predictions 
IZhang & MacFadyenl (|2009l) . 

4.2. Late time calorimetery andt s 

Given the good agreement between the recent wide jet si m- 
ulations IZhang & MacFadven! [20091: IWvgoda et ail 1201 lb nd 
that at late times t > ?nr the shape of the light curve is de- 
term ined by rather than f?o, we also support the conclusion 
from lWygoda et all (1201 ll) that using a Sedov-Taylor approx- 
imation based on Ej to model the light curve at late times 
is a valid approach. For observers approximately on-axis, 
and when errors on the flux of the order of a few are accept- 
ab le, this approach is ev en viable at times t > f,. m , as defined 
by | Wvgoda et alJ (1201 ll) (e.g. for the analysis in iFrail et all 
2000). But this does not mean that the jet is approximately 
spherical b y this time, and neither does it confirm the as- 
sumption in lLivio & Waxm an (2000) of fast deceleration over 
a neglible increase in distance. 

The latter is illustrated in Fig [5] It shows the on-axis blast 
wave radius both for the narrow and wide jet simulation. 
According to Wygo daet al.l d201 ll) . f Siffi is given by eq. [3j 
where the increase in radius after tg is neglected. As a result, 
f 5 .0 « 148 days for the narrow (6a = 0.05 rad) jet. However, 
if t s < t^R the increase in radius can not be neglected com- 
pared to the increase in time since to, because 7 > 1 implies 
that R(t s ) = R s will be comparable to ct s . From Fig [5] we 
find r; e = t s -R' s = 295.86-7.64455 ■ 10 17 cm/c = 0.71 days, 
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FIG. 5. — On-axis blast wave radius as a function of lab frame time for 
narrow and wide jets. The vertical lines denote from left to right: tg (onset of 
spreading), t s ~ 2tg (spreading originally expected to be completed), t(-y = 2) 
(IBM approximation ceases to be applicable, deceleration expected to start to 
become apparent even for spherical jets), tf/n (jet in non-relativistic regime). 

where primed quantities are based on the simulation results 
rather than analytically estimated. For comparison, if there 
were no additional deceleration due to sideways spreading at 
all, and the blast wave had kept on decelerating according to 
7 = j e {t/teT i/2 (BM), the resulting t'/ @ = 0.37 days. Clearly, 
0.71 days lies closer to 0.31 days than to 148 days. For the 
wide blast wave t ' sm = 1 5.8 days, while t s> Q = 372 days. 

Wyg oda et al.l (1201 11) show that beyond f,®, the differ- 
ence between the simulation light curve and a simple non- 
relativistic model is no more than a factor of a few. And 
indeed, although the above demonstrates that it is not valid 
to interpret this as confirming the (modified) Rhoads model 
of quick relativistic spreading and deceleration, t s ,® remains 
useful as a mle-of-thumb to indicate the point when the non- 
relativistic model starts to become appropiate for calorime- 
try estimates (but depending on the amount of error that is 
deemed acceptable and only for observers close to the jet axis, 
see below). Mathematically compensating the lack of quick 
spreading by neglecting the increasing radius has the effect 
that t/ Q is put far into the transrelativistic regime (even for 
spherical blast waves) and as such it has the advantage of 
marking the turnover into the semi-spherical non-relativistic 
phase, while improving upon t^R by retaining a dependency 
on 6j. Insofar as a clear interpretation of f Sl ® is possible, it 
is that f, ® reveals the effects of E: and logarithmic spreading 
throughout the transrelativistic regime. In a practical sense, 
studies that dep end on t s m remain valid for rough calori metry 
estimates (e.g. iFrail et alfeoOOtlShivvers & Bergerll201 lb . 

For more accurate calorimetry and in order to model sig- 
nificant late time features such as the rise of the counterjet, 
detailed jet simulations remain indispensable. A straightfor- 
ward illustration of nonsphericity is given in Fig. |6] which 
shows radio light curves at t > t s for observers at high angles. 
For such orphan afterglows, a calorimetry calculation based 
on the assumption of (semi-) sphericity beyond t s is severely 
limited given that the light curves can still differ between ob- 
server angles by more than an order of magnitude, even if the 
rapid variations may render orphan afterglows easier to detect 
by surveys. For observers slightly off-axis the issue is less 
problematic, as shown by Fig. [6] and again, in general, one 
can assume that as long as the prompt emission is observed, 
radio calorimetry will yield approximately correct results. 
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FIG. 6. — Late time radio light curves observed at 1.43 GHz, for a nar- 
row jet with 6»o = 0.05 rad and Ej = 1.25 X 10 50 erg. Synchrotron self- 
absorption is included but hardly affects the light curves. Light curves are 
plotted for various observer angles, including o b s > #o (i- e - an orphan after- 
glow). The vertical dash-dotted grey curve denotes t s , the point at which the 
jet was originally expected to become quasi-spherical. The calculation uses 
z = 0-5, luminosity distance = 8.7 X 10 27 cm, p = 2.5 and the method from 
Ivan Eerten & MacFadven| <20Tlh . 

It should be noted however, that this does not extend to any 
modeling that includes a measurement of the slope of the late 
time radio light curve. Even without effects such as the rise 
of counterjet (that even leads to a temporary increase in flux 
for a narrow jet), the light curve slope has not yet settled into 
the ST regime. This implies that an attempt to determine cir- 
cumburst density parameter k (for a profile p oc r~ k ) and/or 
the synchrotron slope p from t he light curve slope will stil l 
be bound by the late 5?/vr value dZhang & MacFadven 2009), 
rather than f ,.®. Similarly, a flattening of the light curve be- 
yond f s .® does not automatically imply the rise of a flux con- 
tribution frorn_tJiehost_galaxy (e.g. as assumed for GRB 
98Q703. lBerger et al.ll200lh . 

5. CONCLUSIONS 

We draw the following conclusions with regard to the ques- 
tions posed in sectionQ] 

1 . There is no regime of exponential sideways expansion 
that noticably advances the onset of sphericity, not even 
for narrow jets. The front of the blast wave is far from 
homogeneous, which is a key requirement for the ap- 
proximation by Rhoads (1999) to hold. Even as mate- 
rial at the edges of the jet moves sideways, the bulk of 
the material remains unaffected for a long time. 

2. Sideways expansion of regions of fixed jet energy frac- 
tion is logarithmic in the entire region after the onset 
of expansion (different for different energy fractions) 
and ?nr- The jet remains highly unspherical at ?nr- The 
full transition to sph ericity is a very slow process taking 
~ 5?nr to complete (Zhang & Ma cFadvenl2009t) . 

3. The post jet break light curve slope is determined by 
both spreading and jet edges becoming visible. Com- 
pletely ignoring jet spreading at this stage will notice- 
ably underestimate the steepening of the jet. 

4. Although its impact is slightly decreased by jet spread- 
ing, the observer angle remains important, even for nar- 
row jets and off-axis observation of the jet will result in 
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a light curve where the break time can be significantly 
postponed and even buried in the noise. 



5. Calorimetry based on applying a spherical nonrelativis- 
tic jet model to radio afterglows will remain relevant 
and approximately correct beyond t s when the ob- 
server is approximately on-axis and errors of a few in 
flux level are considered acceptable. However, this does 
not imply sphericity of the outflow and therefore does 
not translate to high observer angles relevant to orphan 
afterglows. Also, t Sj Q represents the effects of Ej and 
long-term spreading during the transrelativistic phase, 
rather than quick early time spreading and decelera- 
tion (because the original derivation of t s ,® based on 
the latter underestimates the increase in radius during 
the spreading phase). For more accurate calorimetry 
and in order to model significant late time features such 
as the rise of the counter] et and features that depend 
on the light curve slope, detailed jet simulations remain 



indispensable. 

6. The magnetar model is challenged less by recent obser- 
vations once allowance is made for the fact that the ob- 
server angle is not known when jet opening angle and 
energy are inferred from the jet break time. The in- 
ferred energy decreases by a factor 2-4. 
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